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ABSTRACT 


The  frequency  dependence  of  the  electrolytic  conductivity 
was  studied  for  solutions  of  varying  concentrations  of  NaCl , 
MgS04,  KC1,  and  KBr.   An  experimental  test  fixture  was 
designed  and  an  equivalent  electrical  model  of  the  test 
system  developed.   A  theoretical  model  of  the  conductivity 
that  accounts  for  charge  carrier  inert: a  is  proposed. 
Measured  values  of  the  impedance  at  various  frequencies 
were  used  to  generate  test  system  model  parameters,  and 
subsequently  identify  sample  response.   Interpretation  of 
the  sample  response  using  the  conductivity  model  allowed 
determination  of  the  conductivity,  which  is  presented  in 
the  form  of  the  D.C.  value  or  real  part,  Ko  and  the  non- 
dielectric  time  constant,  Tc.   The  conductivity  of  the 
solutions  decreased  with  increasing  frequency  and  the 
initial  measurements  of  Tc  were  of  the  order  of  nanoseconds. 
Variations  in  K  with  salinity  were  in  agreement  with  the 
empirical  formula  of  Walden. 
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I.   INTRODUCTION 

Extensive  research  has  been  directed  at  examining  the 
propagation  of  electromagnetic  radiation  in  all  environments. 
0£  fundamental  concern  to  the  Navy  has  been  propagation 
through  its  environment;  the  atmosphere,  the  air-ocean 
interface  and  the  ocean  itself.   Investigations  concerning 
ocean  electromagnetics  require  an  understanding  of  seawater 
conductivity,  or  more  fundamentally,  electrolytic  conductivity 
Historically,  electrolytic  conductivity  research  has  been 
divided  into  two  areas;  low  frequency  research  using 
frequencies  below  10kHz,  and  high  frequency  experiments 
utilizing  frequencies  above  lOMhz  which  probe  the  dielectric 
nature  of  the  solutions. 

In  this  thesis,  the  frequency  dependence  of  the  electro- 
lytic conductivity  is  observed  between  10kHz  and  10MHz  by 
measuring  its  impedance.   The  salt  solutions  examined  are 
the  dominant  contributors  to  seawater  conductivity. 
Sections  II  and  III  present  current  theory  and  model 
development.   Sections  IV  and  V  describe  the  experiment 
and  present  an  analysis  of  the  results.   Section  VI 
concludes  with  a  summary  and  some  possible  areas  for  future 
investigation. 


II.   BACKGROUND 

An  extensive  literature  search  was  carried  out  utilizing 
general  chemistry  reference  material,  textbooks,  and 
computer  methods.   The  Chemistry,  Ocean  Sciences,  Geophysics 
and  Electrical  Engineering  data  banks  of  the  DIALOG  (Lockheed 
Data  Base)  computer  information  system  contain  information 
dating  from  approximately  1967.   Systematic  pursuit  of  the 
reference  materials  and  their  associated  bibliographies/ 
references  consistently  lead  to  the  Debye-Falkenhagen  model 
of  conductivity  proposed  in  1928  [Ref.  1].   Harned  and  Owen 
[Ref.  2]  (1963)  along  with  Condon  and  Odishaw  [Ref.  3]  (1963) 
state  Falkenhagen' s  theory  and  its  historical  verification. 
Smedley  [Ref.  4]  (1980)  describes  improvements  to  the  Debye- 
Falkenhagen  theory  that  occurred  in  the  1970' s,  but  these 
newer  theories  apply  to  low  concentration  solutions  in  the 
low  frequency  regime  (less  than  10  kHz)  only.   Thus,  the 
most  recent  model  that  describes  the  excitation  frequency 
dependence  of  conductivity  was  published  in  1928. 

The  model  for  Falkenhagen' s  conductivity  theory  is  that 
of  a  hard  sphere  ion  under  the  influence  of  an  applied  electric 
field  drifting  in  a  viscous  and  permeable  medium  [Ref.  5]. 
The  reduction  in  mobility  as  concentration  is  increased  is 
due  to  coulombic  interactions  between  ions,  the  predominant 
effects  being  relaxation  and  electrophoresis. 
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An  ion  of  charge  Q  is  surrounded  by  ions  whose  net  charge 
is  -q.   When  an  external  electric  field  is  applied,  the 
central  ion  will  be  attracted  to  the  electrode.   Consequently, 
the  previously  spherically  symmetric  field  around  the 
central  ion  becomes  asymmetric,  for  it  cannot  "relax"  fast 
enough  to  follow  the  motion  of  the  central  ion  toward  the 
electrode.   This  creates  a  small  force,  the  relaxation  force, 
which  inhibits  the  mobility  of  the  central  ion.   When  an 
ion  moves  in  an  electrolytic  solution,  it  tends  to  drag 
the  local  solvent  molecules  with  it.   Since  anions  and 
cations  move  in  opposite  directions,  each  ion  is  effectively 
moving  against  a  stream  of  solvent  molecules.   The  subsequent 
reduction  in  ion  mobility  by  this  effect,  the  electrophoretic 
effect,  can  be  attributed  to  an  electrophoretic  force. 

Falkenhagen  has  included  the  force  due  to  the  applied 
field,  the  Stoke 's  law  hydrodynamic  force  for  a  hard  sphere 
in  a  viscous  continuum,  the  electrophoretic  force  and  the 
relaxation  force  in  developing  an  equation  of  motion  for  the 
ion.   He  does  not  include  any  inertia  term  and  stated  "...it 
is  permissible  to  neglect  the  forces  due  to  dynamical 
reactions  in  comparison  to  the  viscosity  forces."   [Ref.  6]. 
The  forces  due  to  dynamical  reactions  represent  the  inertia 
of  the  ion  in  the  solution,  i.e.  its  inherent  resistance 
to  a  change  in  velocity.   Falkenhagen  qualified  his  statement 
by  comparing  the  two  forces  and  showing  that  the  viscosity 
force  was  larger  than  the  inertial  force.   Two  assumptions 
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critical  to  his  argument  concern  the  charge  carriers, 
specifically  the  size  (radius)  and  mass.   In  view  of  present 
uncertainty  as  to  the  methods  and  mechanisms  of  charge  transport, 
the  validity  of  Falkenhagen' s  assumptions  are  suspect. 
The  absence  of  this  inertia  term  lead  to  a  real  expression 
for  the  conductivity.   Falkenhagen' s  relaxation  force  is 
inversely  proportional  to  the  frequency  of  the  applied 
electric  field,  and  decreased  with  increased  frequency. 
Consequently,  he  stated  that  conductivity  will  increase  with 
increasing  frequency  for  the  ion  is  more  mobile.   Experimental 
data  supporting  this  theory  was  collected  by  Sack  [Ref.  7] 
and  other  investigators  and  summarized  by  Geest  [Ref.  8], 
Careful  examination  of  these  experiments,  all  performed  in 
the  late  1920's,  revealed  that  they  were  relative  measurements 
obtained  by  recording  the  difference  in  response  of  both 
the  test  cell  and  solution  at  various  frequencies  using  a 
bridge  network.   Additionally,  the  measurements  were  subject 
to  sizeable  experimental  uncertainties  (of  the  same  magnitude 
as  the  measured  parameter)  due  to  the  sensitive  nature  of 
the  measurement  and  the  available  technology.   Thus,  the 
results  from  these  experiments  required  careful  interpretation 
and  were  based  heavily  on  existing  theory. 

Retaining  the  inertia  term  in  the  equation  of  motion 
results  in  a  complex  form  for  the  conductivity,  and  is 
developed  in  section  III.   As  the  frequency  is  increased, 
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the  imaginary  part  of  the  conductivity  becomes  smaller.   The 
net  effect  is  a  decrease  in  conductivity  with  frequency. 
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III.   THEORY 

This  section  presents  the  theoretical  foundation  upon 
which  the  experiment  was  conducted.   Part  A  presents  the 
fundamental  development  starting  from  Maxwell's  equations. 
Part  B  develops  the  expression  for  the  conductivity  and 
Part  C  describes  the  specific  application  to  this  experiment. 

A.   FUNDAMENTAL  BACKGROUND 

Maxwell's  equations  for  time  harmonic  fields  can  be 
expressed  as : 

7  x  B  =  li  J  +  ujoiD   ;  V  x  E   =  -joiB 

V  x  D  =  p   ;  V-B  =  0 

the  constituative  relationships  are: 

J  =  KE   ;   B  =  yH 

D  =  eE  =  £0E~  +  P   =  £0E   +  e0x  E  =  e0(l    +  x)  E   where, 
K  is  the  electrical  conductivity,  X  =  x'  -  jx"  is  the  complex 
susceptibility  and  e  =  e   (1  +  X)  is  the  complex  permativity. 
After  substitution,  Maxwell's  equations  reduce  to: 

V  B   =  yo  (K  +  joje)  E  =  \io   j  e  E 

V  E   =  -jeB  ;VB  =  0 
V-D  =p   ;  e  -  e  (1  +  S — ) 

JJ   =eE;   p  =  £V-E=   eV-E+   £-  V  •  E 

J  CO 
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where  the  effective  quantities  (denoted  by  v)  have  been 
introduced  for  computational  ease.   This  formulation  is 
typically  used  in  problems  involving  conducting  dielectrics. 
To  more  clearly  distinguish  between  dielectric  and  conductor 
properties,  several  specific  cases  are  examined. 
1.   Case  1 

Consider  a  conducting,  lossy  dielectric  medium 
between  two  parallel  plates.   Using  the  effective  notation, 
the  impedance  of  this  device  is : 

z   =  - — —  where 
J^c 


c  = 


sE  •  ds    e  A       A  =  plate  surface  area 
E  '  d£  I   =  plate  separation 


and      e      =  e  (1  +  - ) 

Therefore  the  impedance  (in  conventional  notation)  is 

1  1 


z  = 


joiey         jo)  A(e-iJL) 


This  result  could  be  generated  by  a  more  conventional  approach 
The  medium  can  be  characterized  by  a  capacitor  in  parallel 
with  a  resistor.   The  impedance  is: 

where   c   =  -r-        and 


so  that 


j  qjc+G  I 

r  KA 

G   "  1T 


1 


•ojc  A   KA       7X   T~K~, 
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.  A,  .K. 


which  is  identical  to  the  earlier  result.   Incorporating 
the  constituative  equations,  which  are: 

e  =  e0  (1  +X)  =  e'-  je"   and 
X  =  x    -  jx" 


the  impedance  of  the  device  can  be  expressed  as: 


z  = 


1 


JG>f    £V+£(CO£0X"+     K) 


(1) 


2.   Case  II 


Consider  a  finite  conductor  between  two  parallel 
plates.   The  impedance  is  purely  resistive  and  of  the  form 


z  =  R  = 


AK 


C2-)     £  ■  length  of  conductor 

A  =  conductor  cross -sectional 
area 

K  =  electrical  conductivity 


3.   Case  III 


Consider  a  non-conducting  lossy  dielective  between 
two  parallel  plates.   The  impedance  is: 
1  1 


z  = 


JO)C 


j*>  t(£'-  je") 


(3) 


where 


=  JsE-ds    A    (    ,     ,nA 
J  E-d£     *  * 
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B.   CONDUCTIVITY  MODEL 

A  standard  definition  for  conductivity  is: 

K  =  p+  u+   +  p-  u-  (4) 

where  p+,-  is  the  charge  density  of  the  +,-  ion  and 

u+ , -  is  the  charge  mobility  of  the  +,-  ion.   The  mobility 

of  the  ion  can  be  obtained  following  Jackson  [Ref.  10]: 

<v>  rm 

m  ^~  +  mb<v>  =  qE  (6) 

where  <v>  is  the  mean  velocity  of  the  ion,  m  is  the  mass, 
b  is  a  damping  constant  that  reflects  the  change  in  ion 
velocity  due  to  collisions,  q  is  the  charge  of  the  ion, 
and  E  is  the  applied  electric  field.   Only  one  dimensional 
motion  will  be  considered.   Let  <v>  and  E  have  the  time 
harmonic  form  exp(joot).   Then,  solving  equations  (5)  and  (6) 

we  obtain 

Uo  ,      rT    q  (7) 

u   =   -,  r  •  f — nrv     where   Uo  =  ^^-  ^    J 

I  +  jCoj/dT  mo 

Substituting  equation  (7)  into  equation  (4) ,  the 
expression  for  conductivity  becomes: 


K  =   p+Uo+     +   p-Uo- 

1  +  3u+  1  +  Jk- 


(3) 
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Considering  the  simple  salts  of  interest,  such  as  NaCl , 
we  can  assume  p+  =  p-  =  p.   Additionally,  given  equal  anion- 
cation  charge  magnitudes,  the  damping  forces  are  assumed 
to  be  similar.   Therefore,  b+  =  b-  =  b.   Applying  these  two 
simplifications,  the  expression  for  the  conductivity,  from 
Becker  [Ref.  11],  becomes: 

K  =  [p(Uo+  +  Uo-) ]   =    Ko 

,      .  CO  ,      .  CO 

1  +  Jb  1  +  ^ 

where    Ko  =  p(Uo  +  Uo-) 

C.   APPLICATION 

The  experiment  dealt  with  simple  salt  solutions  which 
could  be  considered  as  conducting,  lossy,  dielectric  media. 
The  test  fixture,  described  in  Section  III,  can  be  modeled 
as  two  parallel  plates.   Therefore,  the  impedance  of  the 
test  circuit  can  be  interpreted  as  that  derived  in  case  I 
using  equations  (2)  and  (3).   Incorporating  the  expression 
for  the  complex  conductivity,  the  impedance  of  a  conducting, 
lossy,  dielectric  medium  between  paralle  plates  becomes: 

-    f  .coAe '    A    .  ,    AK  ' l 


_  (  .  cuAe ,    A    ,  ,    A    Ko        ^  " l 

-  ()—   +  i  W£     +  i  (i  ;  jco/bj    } 


=  Cjcoe"   +  jcoe'   +    l  .      — ) 

AKo   +   ^AKob 
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i 


1       -1 

=    (icue11      +      jcoC      +      5   r    7    T ) 
w  J  R  +    jcoL/ 


where  «  « 

R  =      AKo"     and   L   =      F 


Since  this  research  considered  only  aqueous  solutions, 
and  utilized  frequencies  below  10  MHz,  several  simplifications 
can  be  made.   From  Kittle  [Ref.  12],  the  dielectric  constant 
of  water  at  room  temperature  is  given  by: 

+   4rraN     ~   8Q 

l  +  Ot)2 

where  a  is  the  static  orientational  polarizability ,  t  the 
dielectric  relaxation  time  =  10"11  ,  and  N  the  number  of 
spherically  shaped  molecules.   For  co  =  1C8  ,  we  have  (^t)  = 
10"3  so  that: 

e'  =   1  +  lVao-TT2  s   l  +  47TaN  s  80 

and  is  independent  of  frequency.    Thus,   4  tt  aN  =  79. 

Also  from  Kittel  [Ref.  13],  the  complex  dielectric  constant 

is : 

'  ■    t  i  i         4tt    aN  .  4  it  aN    t 

e  Z         '      3e  -      1     +      1     +      (at)  2      "      Jl     +      (0)t)2 


If  we    compare    the    real    and    imaginary  parts: 


£ 


4  7T  aN(cot) ~  79  (cot) 


£ 


1    +    (art)2    +    47raN  1    +    (cot)2    +    79 
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79  (<ut) 

80  +    (cot)2 


Now,    for  a)  _<10    ,    and      t   =    10""s;    the    ratio   becomes 


e"         79(1Q-U)  ~      79(10'3) 

e'         80   +    10"22co2  80   +    CIO"6) 


=      10"3 

As   decreases,  e"  gets  smaller  much  faster  than  e',  and 
the  difference  between  them  grows.   Therefore,  we  can 
neglect  the  e"  term  in  the  impedance  equation. 
The  impedance  becomes: 


Z  =   UOJC  +  -pr— — : — r) 


This  form  of  the  impedance  corresponds  to  the  circuit  model 
shown  in  Figure  1.   The  capacitor  C  reflects  the  effects 
of  the  solution  in  the  conductivity  cell.   The  resistor  R 
and  inductor  L  reflect  the  effect  of  the  specific  salt 
chosen  as  solute.   The  experiment  consisted  of  measuring 
the  value  of  C  using  distilled  water  and  measuring  the 
total  impedance  of  the  cell  filled  with  sample.   From  this 
data,  the  values  of  R  and  L  were  calculated,  and  subsequently 
the  value  for  b  (b  =  R/L) . 

It  is  important  to  note  that  the  e "  term  which  was 
neglected  in  the  impedance  equation  would  correspond  to 
a  resistive  effect  in  light  of  the  proposed  circuit  model. 
It  affects  the  real  part,  not  the  imaginary  part  of  the 
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Figure  1. 


Electrical  Equivalent  Circuit 
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impedance.   Thus,  the  sign  of  the  phase  angle  (negative 
implying  capacitive  behavior;  positive  implying  inductive 
behavior)  would  be  determined  by  the  other  circuit  elements 
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IV.   EXPERIMENT 

A.   APPARATUS 

1 .   Measurement  System 

Impedance  measurements  were  made  with  the  Hewlett- 
Packard  (HP)  Multi-Frequency   LCR  Meter  (type  4275A, 
HP  No.  2045J0]046)  on  loan  from  the  Applied  Physics  Laboratory, 
The  Johns  Hopkins  University.   This  instrument  was  calibrated 
by  Hewlett-Packard  18  November  1982  using  calibration  standards 
traceable  to  the  National  Bureau  of  Standards  to  the  extent 
allowed  by  the  Bureaus'  calibration  facilities  [Ref.  14]. 
Control  of  the  LCR  meter  was  through  a  Hewlett-Packard 
85  Personal  Computer  (HP  No.  2139A4139A)  via  the  Hewlett- 
Packard  Interface  Bus . 

The  LCR  meter  is  a  microprocessor  based  impedance 
measuring  instrument,  [Ref.  15],  which  measures  the  vector 
impedance  (or  admittance)  of  the  unknown  sample  to  be 
tested.   Ten  test  frequencies  were  available  from  the  LCR 
meter:   10kHz,  20kHz,  40kHz,  100kHz,  200kHz,  400kHz,  1MHz, 
2MHz,  4MHz,  and  10MHz. 

Connection  of  an  unknown  sample  was  as  shown  in 
Figure  2.   A  four  terminal  (HPOT,  HCUR,  LPOT,  and  LCUR  on 
Figure  2)  network  was  used  to  connect  the  LCR  meter  to 
the  device  under  test  (DUT) .   This  terminal  architecture 
limited  the  effects  of  mutual  inductance,  interference  of 
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the  measurement  signals,  and  unwanted  residual  factors  in 
the  connections  which  are  normally  encountered  in  high 
frequency  measurements.   The  measurement  current  utilized 
the  outer  shield  conductor  for  a  return  path.   Since  the 
same  current  flowed  through  inner  and  outer  conductors , 
but  in  opposite  directions,  no  net  inductive  magnetic  field 
was  formed  which  ensured  minimal  error  contribution  by  the 
test  leads  or  fixture  to  the  measurement. 

Dependence  of  the  measurement  on  the  test  fixture 
was  also  minimized  by  the  use  of  the  LCR  meter  zero  offset 
adjustment  (ZOA)  [Ref.  15].   The  inductive  and/or  capacative 
nature  of  the  complete  test  fixture  was  measured  at  each 
frequency.   This  was  done  in  two  steps,  first  by  measuring 
the  capacitance  and  conductance  of  the  fixture  in  an  open 
circuit  state  (e.g.  empty,  dry  conductivity  cell).   Then 
the  impedance  and  resistance  of  the  fixture  in  the  short 
circuit  state  (e.g.  conductivity  cell  filled  with  mercury) 
were  measured  at  each  frequency.   The  meter  retained  these 
values  and  automatically  performed  optimum  compensation  on 
subsequent  measurements  to  remove  test  fixture  response. 

The  HP- 85  computer  was  used  to  control  the  LCR  meter 
via  the  interface  bus.   Several  BASIC  language  programs 
were  written  that  fully  controlled  the  measurements  taken 
on  a  given  sample.   This  procedure  ensured  that  the  LCR  meter 
setup  (i.e.  measurement  parameter,  test  signal  level, 
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frequency,  etc.)  was  identical  for  each  test  data  cycle.   It 
also  enabled  a  large  number  of  measurements  of  a  given 
parameter  to  be  accomplished  in  a  short  amount  of  time. 
2.   Test  Cells 

Two  test  cells  were  used  in  this  experiment.   Cell  #1 
was  a  conventional  conductivity  cell  as  sketched  in  Figure  3. a 
The  borosicilate  glass  cell  body  held  approximately  30ml  and 
had  2  electrodes  coated  with  platinum  black.   The  electrode 
leads  penetrated  the  cell  wall  through  lime  glass  supports 
and  a  brass  capped  connection  point.   The  cell  constant 
is  defined  as  the  ratio  of  the  separation  distance  to  the 
surface  area  of  the  electrodes.   The  cell  #1  cell  constant 
of  1.19  1/cm  was  determined  using  the  procedures  of 
Reference  16.   A  solution  consisting  of  .7466g  KC1  in  1KG 
of  solution  (KC1  plus  H20)  was  the  standard.   The  resistance 
was  measured  using  a  LKB-PRODUKTER  conductivity  bridge 
(type  LKB  3216B)  calibrated  18  March  1983. 

Test  cell  #2  is  depicted  in  Figure  3.b.   The  cell 
was  made  of  5/16  inch  ID  thickwall  TYGON  tubing  and  two 
nylon  "tee"  connectors  which  provided  fill  and  drain  ports. 
The  electrodes  were  made  of  5/16  inch  OD  solid  carbon  rods 
machined  such  that  the  electrode  face  was  flat  and  perpen- 
dicular to  the  axis  of  the  rod.   The  electrodes  were 
connected  to  the  instrument  test  leads  by  means  of  machined 
brass  clamps  as  shown  in  Figure  4.   The  cell  constant  of 
27.38  1/cm  for  cell  #2  was  determined  by  comparison  to 
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Figure  3.   Experimental  Conductivity  Test  Cells 
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Figure  4.   Electrode  Connection  Clamp 
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cell  cell  #1;  the  value  of  conductivity  for  a  specific 
solution  must  be  the  same  for  both  cells.   This  method  was 
used  because  cell  #2,  by  virtue  of  its  design,  cannot  be 
thermostated  as  accurately  as  cell  #1. 

Each  cell  was  supported  and  partially  enclosed  in 
styrofoam.   This  support  arrangement  ensured  no  additional 
electrical  or  magnetic  pertubation  of  the  measurement 
and  minimized  any  thermal  fluctuations.   Temperature  measure- 
ments were  made  with  a  WEKSLER  (type  1509)  immersion  thermo- 
meter.  The  experiment  was  performed  in  a  relatively  static 
thermal  environment.   Early  measurements  of  a  variety  of 
solutions  showed  no  appreciable  temperature  variation 
during  a  measurement  cycle.   A  constant  temperature  bath 
was  not  utilized  because  of  the  static  thermal  environment 
and  the  inductive  effects  observed  due  to  bath  operation/design 

B.   EXPERIMENTAL  APPROACH 
1 .   Experimental  Model 

Determination  of  the  solution  parameters  of  interest 
required  that  any  effects  due  to  the  cell  (test  fixture) 
be  understood  and  eliminated  if  possible.   To  this  end 
the  problem  was  divided  into  two  parts:   the  physical  test 
fixture,  and  the  equivalent  electrical  circuit  simulating 
the  test  solution  and  its  interaction  with  the  test  fixture. 

a.   Physical  Test  Fixture 

The  test  fixture  consisted  of  the  conductivity 
cell  and  the  leads  connecting  it  to  the  measuring  instrument. 
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Since  the  conductivity  cell  consisted  of  two  flat  electrode 
plates  with  a  dielectric  between  them  (primarily  water) 
it  was  capacitive  by  nature.   The  test  leads  were  approximately 
6  inches  each  of  RG-58C/U  coaxial  cable,  and  had  a  small 
resistance  and  inductance  associated  with  them.   Additionally 
there  were  several  coaxial  connector  joints  and  solder 
joints,  each  with  some  residual  effect.   These  accumulated 
effects  were  quite  complex,  but  were  essentially  negated 
using  the  zero  offset  adjustment  (ZOA)  feature  of  the  LCR 
meter.   The  ZOA  was  performed  sequentially  in  two  steps. 
The  test  fixture  was  assembled  with  the  conductivity  cell 
dry  for  the  open  circuit  portion.   Initiated  by  the  "open" 
button,  the  meter  automatically  measured  the  capacitance 
and  conductance  at  each  test  frequency.   The  cell  was 
carefully  filled  with  mercury  and  the  short  circuit  portion 
of  the  ZOA  initiated  using  the  "short"  button.   The 
instrument  automatically  measured  the  inductance  and  resistance 
at  each  test  frequency.   The  values  were  retained  by  the 
LCR  meter  and  subsequent  measurements  were  compensated  to 
remove  test  fixture  response.   To  identify  and  verify  the 
remaining  background  of  the  test  fixture,  the  inductance, 
resistance  and  vector  impedance  of  the  mercury  were  measured 
and  recorded  (see  Appendix  A) . 

Unfortunately,  cell  #1  did  not  have  a  drain 
port.   It  had  to  be  physically  disconnected  to  be  drained 
and  rinsed  for  each  new  test  sample.   This  reconnection 
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altered  the  test  fixture  slightly,  and  could  have  introduced 
errors  not  accounted  for  in  the  initial  ZOA.   To  identify 
these  errors,  the  inductance,  resistance,  and  vector 
impedance  of  mercury  was  measured  after  a  reconnection  and 
compared  to  readings  after  ZOA.   As  shown  in  Appendix  A, 
post  reconnection  values  were  slightly  greater  than  post 
ZOA  values  and  were  subsequently  used  to  determine  measurement 
accuracy.   The  open  circuit  portion  of  the  ZOA  was  performed 
prior  to  each  measurement  run  to  partially  compensate  for 
the  reconnection  change. 

b.   Test  Solution  Equivalent  Circuit 

The  equivalent  electrical  circuit  of  the  test 
solution  shown  in  Figure  5  was  derived  using  the  following 
considerations.   The  capacitor  C  reflected  the  general 
character  of  the  conductivity  cell;  a  parallel  plate  capacitor 
with  a  dielectric  material  between  the  plates.   The 
dielectric  material  was  the  test  solution,  which  was  pre- 
dominantly water  even  at  the  higher  concentrations.   From 
Hasted  [Ref.  17],  the  dielectric  constant  does  depend  on 
salinity,  but  this  is  a  small  effect  and  was  subsequently 
ignored.   Therefore,  the  value  of  C  was  determined  reasonably 
accurately  from  the  measured  response  of  pure  water  as  test 
solution. 

The  rest  of  the  circuit  in  parallel  with  C 
represented  the  response  of  the  solute   and  its  interaction 
with  the  electrodes.   The  electrodes  were  interface  surfaces 
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Figure  5.   Equivalent  Electrical  Circuit 


for  the  transition  between  electronic  conduction  and  ionic 
conduction,  and  on  a  microscopic  scale  extremely  complex. 
However,  on  a  macroscopic  scale,  this  transition  effect 
was  more  simply  modeled  as  a  capacitor  in  parallel  with  a 
resistor.   The  capacitor  reflected  the  dielectric  layer 
which  plated  on  the  electordes,  while  the  parallel  leakage 
resistor  reflected  the  imperfections  in  the  layers  and 
subsequent  non- ideal  capacitance.   Since  there  were  two 
electrodes,  each  with  a  different  ion  layer  structure,  each 
electrode  was  considered  separately.   CI,  Rl ,  C2,  and  R2 
were  the  model  parameters  that  described  the  two  electrodes. 
The  remaining  two  elements,  R  and  L,  represented 
the  response  of  the  solute.   As  described  earlier,  R  and  L 
formed  the  expression  for  the  complex  conductivity,  and 
were  the  parameters  of  primary  interest  in  this  research. 
2 .   Measurement  Procedure 

The  pure  water  and  stock  solutions  were  allowed 
several  days  to  thermally  stabilize  in  the  laboratory.   The 
LCR  meter  was  energized  4-6  hours  prior  to  any  data  runs. 
Just  prior  to  a  series  of  measurements,  the  ZOA  was  performed 
as  follows: 

1.  Conductivity  cell  was  drained  and  connected  to  the 
meter. 

2.  Open  circuit  portion  of  the  ZOA  performed. 

3.  Cell  filled  with  mercury  and  ZOA  completed. 
Individual  measurements  were  made  using  the  following  sequence 
of  steps: 
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a.  Cell  was  disconnected,  emptied,  rinsed  with 
pure  water  and  reconnected  to  the  LCR  meter, 

b.  Open  circuit  portion  of  ZOA  performed, 

c.  Cell  filled  with  sample, 

d.  Temperature  of  sample  measured. 

e.  Test  solution  description  and  temperature 
entered  into  the  HP-85  computer, 

f.  The  measurement  program  was  run. 

This  procedure  was  devised  using  cell  #1.   As  cell  #2  had 
a  drain  port,  draining  and  rinsing  the  cell  was  accomplished 
without  disconnection  of  the  cell.   Therefore  cell  #2  runs 
followed  the  same  procedure  except  for  disconnecting  the 
cell  and  performing  the  open  circuit  portion  of  the  ZOA. 

The  measurement  program  provided  the  specific 
instructions  to  the  LCR  meter.   The  program  was  written 
using  the  LCR  meter  and  computer  operating  manuals  [Ref.  15 
and  18]  and  followed  suggested  sample  programs.   Initially, 
the  program  defined  the  instrument  circuit  mode,  test 
signal  level,  measurement  range  and  trigger  source.   A 
repetitive  measurement  sequence  followed.   At  each  test 
frequency  the  impedance  (magnitude  and  phase),  inductance, 
resistance,  voltage,  and  current  were  each  measured  100 
times,  averaged,  and  recorded.   This  constituted  one 
measurement  run  and  was  repeated  for  each  new  sample. 

C.   SAMPLE  PREPARATION 

Solutions  of  magnesium  sulfate,  sodium  chloride, 
potassium  bromide,  and  potassium  chloride  were  prepared 
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from  pure  solid  chemical.   Since  magnesium  sulfate  is 
deliquesent ,  the  procedure  for  preparation  of  solutions  of 
the  proper  salinity  was  somewhat  complex.   Mallinkrodt 
Analytical  Reagent  Anhydrous  magnesium  sulfate  was  placed 
in  a  clean  dry  pyrex  beaker  and  heated  in  an  electric 
furnace  at  130C  for  several  hours  after  which  the  powder 
was  placed  in  a  Scheibler  lime  glass  desiccator  and  allowed 
to  cool  to  room  temperature.   Stock  solutions  ranging  from 
a  salinity  (denoted  S)  of  .1  to  100  were  desired.   Salinity 
is  defined  as:   100  (wt  of  solute)  /  (wt  of  solution). 
Preparation  of  the  specific  solutions  was  done  by  estimating 
the  volume  of  solute  needed,  obtaining  that  amount  from 
the  desiccator  and  rapidly  weighing  the  solute  plus  polystyrene 
balance  pan.   All  weighings  were  done  on  a  Sartorius 
analytical  balance  type  2403  accurate  to  .0001  grams. 
Several  practice  weighings  were  done  in  this  manner,  and 
the  rate  of  water  vapor  absorption  by  the  magnessium 
sulfate  was  estimated.   Using  these  figures,  the  sample 
weights  were  assigned  the  accuracy  of  .001  grams.   Once 
weighed,  the  magnesium  sulfate  was  transfered  to  a  clean 
dry  volmetric  flask.   Pure  water  was  added  to  the  flask 
using  a  50  ml  precision  burette  until  the  solution  volume 
was  250+/-. 12ml  (as  indicated  by  the  calibration  line  on 
flask) .   The  weight  of  solute  and  water  used  were  recorded. 
Mallinckrodt  U.S. P.  grade  sodium  chloride,  Mallinckrodt 
analytical  reagent  grade  potassium  chloride  and  Fisher 
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certified  research  grade  potassium  bromide  were  used  to 
prepare  the  stock  solutions  of  each  specific  reagent.   Solution 
preparation  was  similar  to  that  for  magnesium  sulfate; 
except  that  the  heating  and  rapid  weighing  necessitated  by 
the  hygroscopic  nature  of  the  magnesium  sulfate  was  not 
required  for  the  other  reagents. 

D.   EXPERIMENTAL  RESULTS 

The  experimental  data  is  divided  into  3  sections; 
background  measurements,  cell  #1  data  and  cell  #2  data. 

Prior  to  the  data  runs,  the  ZOA  was  performed  as 
described  earlier.   After  the  ZOA,  the  response  of  the 
measurement  system  was  recorded  for  two  specific  samples; 
pure  water  and  mercury.   These  results  are  presented  in 
Appendix  A,  and  were  used  in  the  following  section  to 
analyse  the  data. 

Thirty-two  samples  were  measured  using  cell  #1  and 
the  data  is  presented  in  Appendix  B  in  tabular  form. 
Select  data  sets  characteristic  of  the  rest  are  presented 
in  graphical  form.   Potassium  chloride  was  chosen  as 
representative,  and  the  results  for  salinities  of  1,  25 
and  100  are  shown  in  Figures  6,  7,  and  8  respectively. 
Presented  in  these  figures  are  the  raw  data  points  and 
for  comparison,  the  response  of  the  experimental  model 
(solid  line).   Potassium  bromide,  sodium  chloride  and 
magnesium  sulfate  at  a  salinity  of  25  are  also  shown  in 
Figures  9,  10,  and  11  respectively. 
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Figure  6.   Impedance  vs.  Frequency  for  KC1  S=1.0 
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Figure    7.       Impendance    vs.    Frequency    for   KC1    S=25 
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Figure    8.       Impedance  vs.    Frequency    for   KC1    S=100 
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Figure  10.   Impedance  vs.  Frequency  for  NaCl  S=25 
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Figure  11.   Impedance  vs.  Frequency  for  MgS04  S=25 


42 


Cell  #2  data  is  also  presented,  in  Appendix  C  in  tabular 
form.   Sodium  chloride  and  magnesium  sulfate  data  at 
salinities  of  25  and  100  characterized  the  data  and  are 
presented  in  Figure  12  through  Figure  15  in  graphical 
format  similar  to  cell  #1  data. 
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Figure  12.   Impedance  vs.  Frequency  for  NaCl  Cell  #2 
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Figure    13.       Impedance    vs.    Frequency    for   NaCl    S=100    Cell    #2 
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Figure  14.   Impedance  vs.  Frequency  for  MgS04  S=25  Cell  #2 
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Figure  15.      Impedence    vs.    Frequency    for  MgS04    S=100    Cell    #2 
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V.   ANALYSES 

The  raw  data  from  the  experiment  consisted  of  the 
magnitude  of  the  impedance  (denoted  by  Z)  and  the  phase 
angle,  which  were  the  response  of  the  test  fixture  and 
sample.   A  simple  equivalent  electrical  circuit  has  been 
derived  in  section  IV  that  simulated  the  test  fixture  and 
sample.   With  this  circuit,  the  raw  data  was  used  to 
calculate  the  effect  of  the  sample.   Part  A  describes  how 
the  sample  response  was  calculated.   Part  B  presents  this 
secondary  data,  and  Part  C  discusses  accuracy  and  sources 
of  error. 

A.   SAMPLE  RESPONSE 

The  equivalent  electrical  circuit  has  been  derived  in 
section  III,  and  is  presented  again  in  Figure  16b.   The 
value  of  C  was  determined  using  the  response  of  pure  water, 
the  predominant  constituant  of  the  test  solution.   This 
left  the  six  remaining  parameters  to  be  determined  for 
each  test  solution.   The  method  chosen  to  do  this  was 
to  simplify  the  equivalent  circuit  based  upon  the  frequency 
range,  and  generate  approximate  values  for  the  parameters 
after  which  the  magnitude  of  the  impedance  was  calculated 
and  compared  to  the  actual  data.   The  model  parameters 
were  then  manually  changed  in  an  iterative  process  until 
close  correlation  between  calculated  and  measured  data  was 

achieved. 
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This  process  was  conveniently  accomplished  using  the 
HP-85  computer  and  a  simple  program  shown  in  Appendix  D. 
The  frequency  range  of  10  kHz  to  10  MHz  was  divided  into 
three  regions;  high,  mid  and  low.   Within  each  region, 
the  raw  data  was  fit  to  a  first  degree  polynomial  (i.e. 
f(x)  =  a  +  bx)  using  the  method  of  least  squares  [Ref.  19]. 
The  various  circuit  parameters  were  obtained  from  the 
coefficients  of  the  polynomial.   Capacitor  C  had  a  very 
large  impedance  compared  to  the  rest  of  the  circuit  at  all 
frequencies.   Since  the  effect  of  the  capacitor  C  was 
small,  for  it  was  a  large  parallel  impedance,  it  was  ignored 
in  the  simplified  circuits. 

In  the  high  frequency  region,  4  to  10   MHz,  the  circuit 
model  simplified  to  that  shown  in  Figure  16b.   Capacitors 
CI  and  C2,  although  different,  represented  a  small  impedance 
and  short  circuited  resistors  Rl  and  R2 .   The  simple  series 
RL  circuit  remained.   The  raw  data  for  this  region  was 
fitted  to  the  equation     Z2  =  R2  +  (Leo)  2    and  provided 
values  for  R  and  L. 

The  mid  frequency  region,  200  kHz  to  2  MHz,  has  the 
simplified  circut  model  shown  in  Figure  16c.   The  capacitor 
C  was  larger  than  C2  and  shorted  out  resistor  Rl .   The 
inductor  L  also  represented  a  small  impedance  and  was  ignored 
The  data  was  fit  to  the  equation 

Z2  =  (R  +  R2)2  +  (03C2R2)2  (R2  -  Z2) 
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(this  is  linear  in  coA2) .   Using  the  coefficients  of  the 
polynomial  and  the  previously  calculated  value  of  resistor 
R,  values  of  C2  and  R2  were  computed. 

The  low  frequency  region,  10  kHz  to  100  kHz,  circuit 
model  is  shown  in  Figure  16d.   Capacitor  C2  was  smaller 
than  CI  and  presented  a  large  impedance  in  this  region;* 
it  was  largely  masked  by  the  resistor  R2  and  subsequently 
ignored.   The  inductor  L  represented  a  small  impedance 
and  was  also  ignored.   The  data  was  fit  to  the  equation 

Z2  =  (R  +  Rl  +  R2)2  +  (0)C1R1)2C(R  +  Rl)2  -  Z2) 
Using  the  coefficients  and  previously  calculated  values 
for  R  and  R2 ,  CI  and  Rl  were  calculated. 

The  program  then  computed  the  magnitude  of  the  impedance 
based  on  an  analytical  form  derived  from  the  complete 
equivalent  circuit  model.   Using  the  approximate  parameter 
values  previously  obtained,  impedance  values  were  calculated 
for  each  frequency  and  compared  to  raw  data.   An  iterative 
process  followed  where  new  values  were  entered  manually  for 
each  parameter  to  more  closely  duplicate  raw  data.   The 
final  set  of  six  parameters  represented  the  best  modeling 
of  the  raw  data  and  was  recorded  along  with  the  comparison 
impedance  values  at  each  frequency. 

B.   CALCULATED  RESULTS 

As  developed  in  section  III,  the  conductivity  can  be 
expressed  in  terms  of  the  resistance  R  and  inductance  L. 


50 


PC 


r^V UCIh 


rWi 


H 


oi  pc 


u 


OS 


HH 


CrS  J 


H 


CJ 


CJ 


05 


PC 


h 


CJ 


u 


Figure  16.   Simplified  Model  Circuits 
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The  real  part  of  the  conductivity  at  low  frequencies  is 
Ko  =  k/r,  k  =  cell  constant,  and  is  often  referred  to  in 
the  literature  as  the  D.C.  conductivity.   The  imaginary 
part  of  the  conductivity  is  Ko(o>/b)/(l  +  (w/b)2),  which 
can  be  characterized  by  the  damping  factor  b,  which  equals 
R/L.   The  real  and  imaginary  parts  of  the  conductivity 
are  considered  separately. 

The  damping  constant  b  is  the  reciprocal  of  the 
conductivity  relaxation  time  (non-dielectric) .   This  inter- 
pretation follows  from  the  development  in  Section  III, 
specifically  the  equation: 

md<v>/dt  +  mb<v>   =  qE . 

When   the  applied  E-field  is  turned  off,  the  charge  carriers 
return  to  equilibrium  and  the  mean  velocity  becomes: 

<v>  =  <v(t  =  0)>  exp  (-tb). 

The  conductivity  relaxation  time,  1/b,  is  designated  Tc. 
A  summary  of  Tc  for  all  cell  #1  data  and  select  cell  #2 
data  is  presented  graphically  in  Figure  17.   Tc  data  is 
also  shown  in  tabular  form  in  Table  1. 

Displayed  in  Figure  17,  the  general  trend  of  the  time 
constant  was  to  increase  with  increased  salinity,  regardless 
of  the  reagent.   This  suggested  that  the  effect  of  the 
ion-ion  interactions  were  cumulative;  the  mobility  decreased 
as  the  number  of  ions  increased.   Also  of  note  was  the 
relatively  small  variation  in  time  constant  with  either 
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Figure  17.   Non-Dielectric  Time  Constant  Comparison 
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Table  1.   Time  Constant  Data 


Tc  DATA  (ns) 
CELL  #1 

MpS04        KC1        KBr       NaCl 


.  1 

.516 

.784 

.09 

2.64 

1 

.265 

.899 

,56?6 

1  .00 

4 

1.118 

1.897 

1  .527 

2.  10 

10 

1.972 

2.36 

2.  12 

2.88 

25 

2.351 

2.51 1 

2.31 

3.33 

49 

2.481 

2.416 

2.2^4 

3.4a 

64 

^.04 

1  .028 

2.49 

^.298 

100 

3.14 

?.31 

2.  18 

4.25 

CELL  #? 

S      Mi*S0U  NaCl 

25      3.64  2.98 

100     2.9°  2.83 
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salinity  or  reagent.   A  change  in  salinity  by  100  resulted 
in  a  change  in  time  constant  by  less  than  10.   The  variation 
of  Tc  with  S  showed  definite  structure  unique  for  each 
reagent.   Analysis  of  this  structure  would  require  a  more 
sophisticated  theoretical  model,  and  was  not  pursued  further. 
As  discussed  in  part  C,  accurate  results  for  salinities 
less  than  one  were  fundamentally  more  difficult  to  obtain 
due  to  cell  design.   Consequently,  they  were  at  best  an 
approximation  to  actual  sample  response  and  were  ignored 
in  any  analysis. 

The  parameter  R  related  to  the  real  part  of  the 
conductivity,  specifically  Ko  =  1/r  per  unit  length.   The 
Ko  results  for  this  experiment  are  presented  graphically 
in  Figure  18  and  in  tabular  form  in  Table  2.   Following 
Smedley  [Ref.  20],  Ko  was  expressed  in  units  of  semens/meter 
as  opposed  to  one  of  many  possible  historical  forms. 

Early  models  of  the  concentration  dependance  of  Ko  are 
presented  by  Falkenhagen  [Ref.  21]  and  for  low  concentrations 
are  of  the  form: 

Ko  =  S(A  -  B(S)  ".5) 

where  S  is  the  salinity  and  A  and  B  are  constants.   Accordingly, 
data  points  would  form  straight  lines  of  negative  slope  if 
plotted  as  Ko/S  vs.  (s)  ~.5  as  in  Figure  19.   Calculated 
data  correlated  poorly  to  this  model,  as  expected,  for  the 
low  concentration  assumption  was  no  longer  valid.   Experimental 
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Figure  18.   Conductivity  Comparison 
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Table  2.   Conductivity  Data 


Ko  DATA  (S/m) 


s 

Mf*S04 

KC1 

KBr 

NaCl 

.1 

.408 

.0186 

.0119 

.0268 

1 

.1253 

.  162 

.1116 

.175 

4 

.331 

.609 

.4269 

.635 

10 

.7211 

1.583 

1.029 

1.553 

25 

1.51 

3.825 

2.517 

3.5^6 

49 

2.268 

7.361 

4.869 

6.693 

64 

3.248 

9.542 

6.373 

8.337 

100 

4.139 

15.  115 

10.083 

1^.0^8 
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Figure  19.   Conductivity  Comparison;  Falkenhagen  Model 
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data  more  closely  fit  an  emperical  formula  by  Walden  [Ref.  22] 
of  the  form: 

Ko  =  SA/(1  +  B(S)  ~.5) 

Plotting  S/Ko  vs.  (S)  ~ .  5 ,  data  points  would  form  straight 
lines.   Data  presented  in  Figure  20  approximated  this 
linear  relationship  very  well  at  higher  salinities. 

C.   ACCURACY  AND  SOURCES  OF  ERROR 

This  section  has  been  divided  into  two  parts.   Part  1 
deals  with  the  accuracy  and  errors  associated  with  the 
measurement  process  while  Part  2  discusses  the  errors  in 
the  derived  quantities. 

1 .   Measured  Values 

The  data  consisted  of  impedance  (magnitude  and 
phase)  at  a  specific  frequency  for  a  given  solution  of 
reagent.   Therefore  the  overall  accuracy  was  a  function  of 
the  accuracy  of  each  of  these  three  parameters. 

The  HP-4275A  LCR  meter  was  used  to  measure  the 
impedance.   The  meter  operating  manual  [Ref.  23]  describes 
the  accuracy  of  the  magnitude  of  the  impedance,  designated  Z, 
and  phase  angle  separately.   Z  accuracy  was  a  function  of 
both  the  test  frequency  and  the  measurement  range.   High 
frequency  measurements  were  inherently  more  difficult  to 
make  and  were  subsequently  less  accurate.   Table  3  presents 
a  summary  of  the  impedance  measurement  accuracy.   The 
ten  test  frequencies  were  preset  by  the  LCR  meter  with  an 
accuracy  of  .01%. 
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Figure  20.   Conductivity  Comparison:   Walden  Model 
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Table  3.   Impedance  Accuracy 


PARAMETER  FREQUENCY  RANGE 

VALUE    iOkHz-iOOkHz  200kHz-iMHz   2MHz   4MHz-iOMHz 


Z<i?9 

(OHMS) 

* 

.27. 

TV 

•  -  . . 

.77. 

2.1X 

Z>199 

.4% 

.5% 

1  .  i% 

3.1% 

PHASE 
ANGLE 

.  i 
(DEG) 

.  i 

.3 

.9 

*  .27.    OF  METER  READING  =  POTENTIAL  ERROR 
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Preparation  of  the  samples  involved  two  primary 
sources  of  error;  the  liquid  volume  and  the  sample  weight. 
Volume  measurements  were  made  using  a   50ml  precision  burette 
with  an  accuracy  of  .1ml  for  sample  concentrations  greater 
than  salinity  of  4.   For  lower  concentrations,  liquid 
volume  was  determined  using  volumetric  flasks  accurate  to 
.12ml.   All  sample  weights  were  measured  using  the  Sartorious 
balance  described  earlier,  accurate  to  .OOOlg.   Stated 
salinity  value  accuracy  was  determined  to  be  1%    following 
consideration  of  the  absolute  accuracy,  the  sample  preparation 
procedure  and  potential  changes  in  the  solutions  while 
stored  in  the  laboratory  (i.e.  evaporation/condensation). 

An  additional  source  of  error  not  compensated  for 
by  the  instrument  was  the  skin  effect  inductance  of  the 
electrolytic  solution.   Calculation  of  this  effect  for  an 
electrolytic  solution  and  mercury  are  presented  in 
Appendix  E.   The  ZOA  was  performed  using  mercury,  which 
from  Appendix  E  had  a  skin  effect  inductance  of  about 
2nH.   This  small  inductance  became  part  of  the  compensation 
for  subsequent  readings.   But  the  electrolytic  solutions 
had  a  skin  effect  inductance  of  about  5nH,  and  the 
difference  between  them  (about  3nH)  was  the  skin  effect 
inductance  error.   Although  present,  the  size  of  this  error, 
when  compared  to  the  rest  of  the  equivalent  circuit,  was 
very  small.   The  previously  stated  impedance  accuracy  values 
adequately  account  for  this  additional  error. 
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2.   Calculated  Data 

The  accuracy  of  the  calculated  results  was  a  function 
of  the  accuracy  of  the  raw  data,  the  numerical  manipulation 
of  the  raw  data,  and  the  accuracy  of  the  cell  constant. 

The  determination  of  sample  response  has  been 
described  earlier.   A  first  degree  polynomial  was  fit  to 
the  raw  data  using  the  method  of  least  squares.   From 
the  coefficients  of  the  polynomial  vales  for  all  six 
parameters  in  the  equivalent  circuit  were  calculated.   As 
discussed  in  Hornbeck  [Ref.  24],  the  method  of  least  squares 
inherently  results  in  a  poorly  conditioned  coefficient 
matrix  with  overall  numerical  accuracy  dependent  upon 
the  degree  of  polynomial  used  and  the  number  of  significant 
digits  of  the  computing  device.   Restriction  to  a  first 
degree  polynomial  in  combination  with  the  12  digit  precision 
of  the  HP-85  [Ref.  25]  ensured  a  negligible  calculation 
error  in  the  computed  parameter  values. 

The  final  set  of  parameters  were  determined  via 
an  iterative  process  of  comparing  calculated  impedance 
magnitude,  Z,  to  raw  data  at  each  frequency.   The  accuracy 
of  this  process  was  dependent  upon  the  optimization  criteria 
and  the  sensitivity  of  the  calculated  Z  values  to  changes 
in  any  one  of  the  parameters.   Optimization  consisted  of 
keeping  the  difference  between  calculated  and  measured 
data  to  less  than  1%    (of  the  experimental  value)  at  all 
frequencies.   The  sensitivity  of  the  calculated  Z  values 
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to  variations  in  the  parameters  is  shown  in  Figure  21. 
MgS04  at  S=25  was  chosen  as  a  typical   analysis  case. 
Parameters  were  varied  as  described  from  the  final  parameter 
set,  and  Z  values  recalculated  for  comparison.   The  values 
of  R  and  L  were  directly  related  to  the  test  sample  while 
the  other  four  parameters  were  necessary  to  describe  the 
test  fixture.   Comparison  with  the  other  parameters  showed 
that  the  impedance  associated  with  R  was  much  larger  than 
the  others.   Therefore  the  calculated  Z  values  were  very 
sensitive  to  changes  in  R.   The  value  of  L  was  calculated 
from  data  at  the  higher  frequencies,  for  this  was  the  only 
region  where  the  impedance  associated  with  L  was  significant. 
Consequently,  the  calculated  Z  values  were  sensitive  to 
variations  in  L  at  the  higher  frequencies.   Thus  the  values 
of  R  and  L  served  to  define  the  correlation  between  calculated 
and  observed  data,  and  are  not  present  in  Figure  21. 

Phase  angle  raw  data  was  not  used  in  generating 
the  system  parameters  directly.   For  most  samples,  measured 
phase  angle  values  started  at  near  zero  at  low  frequency 
and  changed  only  a  few  degrees.   The  small  phase  angle 
was  due  to  the  real  part  being  much  larger  than  the 
imaginary  part  of  the  impedance.   The  small  change  in 
phase  angle  reflected  the  small  change  in  impedance  of 
the  inductor  as  frequency  increased.   Thus  the  phase  angle 
was  insensitive  to  changes  in  R  and  L.   Additionally,  the 
measurement  error  associated  with  phase  angle  (Table  3) 
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Figure    21.      Model    Sensitivity    to    Parameter   Variation 
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was  proportionally  much  larger  than  for  impedance  magnitude. 
Therefore  the  small  changes  in  the  phase  angle  had 
incorporated  in  them  a  sizeable  error.   Early  analysis 
efforts  utilizing  phase  angle  comparison  were  discontinued 
due  to  large  variations  in  output  parameters.   The  sensitivity 
and  inherent  error  associated  with  phase  angle  data  relegated 
its  immediate  usefulness  to  that  of  an  indicator  of  data 
trends  only. 

For  samples  of  salinity  less  than  one,  an  additional 
complication  occurred.   When  considering  the  equivalent 
electrical  circuit,  the  parameter  C  was  ignored  in  the 
preliminary  analysis  for  it  was  a  large  impedance  relative 
to  the  rest  of  the  circuit.   For  low  salinity  samples, 
this  relationship  was  not  as  valid;  the  resisitance  R 
of  the  sample  was  much  larger  and  approached  the  impedance 
of  C  within  an  order  of  magnitude.   Thus  the  raw  data 
reflected  the  response  of  C  partially  masking  the  response 
due  to  R  and  L.   As  a  consequence,  the  subsequent  processing 
of  raw  data  to  determine  values  of  R  and  L  was  more 
difficult.   Calculated  impedance  was  much  less  sensitive 
to  variations  in  L  and  R.   Accordingly,  stated  parameter 
values  were  approximations  of  sample  response  and  useful 
as  indicators  of  data  trends  for  this  low  salinity  region. 

The  cell  constant  for  cell  #1  was  determined  to 
be  1.188  +  /-  .0012  1/cm  using  the  data  and  procedure  in 
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Harned  and  Owen  [Ref.  26]  in  conjunction  with  the  LKB 
conductivity  bridge.   The  conductivity  bridge  was  calibrated 
18  March  1983  and  certified  accurate  to  .1%.   Several 
resistance  measurements  were  made  and  the  mean  value  used 
in  computing  the  cell  constant. 


67 


VI.   SUMMARY  AND  RECOMMENDATIONS 

The  central  result  of  this  experiment  was  the  determination 
of  the  electrolytic  conductivity  as  shown  in  Tables  1  and  2. 
The  experiment  was  designed  to  use  direct  measurements 
of  the  impedance  and  included  development  of  a  model  for 
the  electrical  characteristics  of  the  complete  test  fixture 
system  and  a  theoretical  model  for  the  conductivity. 
Isolation  of  the  sample  response  from  the  overall  test 
fixture  response  was  accomplished  using  the  system  model. 
The  conductivity  model  enabled  interpretation  of  the  sample 
response  and  subsequent  calculation  of  the  electrolytic 
conductivity  of  the  sample. 

The  physical  measurement  of  the  electrolytic  conductivity 
was  sensitive  to  the  temperature,  solution  concentration  and 
especially  design  limitations  inherent  in  both  the  test 
fixture  and  measurement  system.   Emphasis  of  the  imaginary 
part  of  the  conductivity  (i.e.  the  inductance  L)  occurred 
near  10  MHz.   The  combination  of  only  preselected  test 
frequencies,  an  upper  frequency  limit  of  10  MHz  on  the  measure- 
ment instrument,  and  a  frequency  dependent  absolute  measurement 
error  that  dominated  the  overall  accuracy  around  10  MHz 
resulted  in  relatively  few  data  points  for  analysis  and 
made  accurate  determination  of  the  inductance  L  more  difficult. 
Additionally,  the  measured  impedance  of  low  concentration 
samples  was  dominated  by  the  capacitive  character  of  the 
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test  cell.   This  largely  masked  the  response  of  the  sample 
and  further  complicated  determination  of  the  sample 
paramters. 

The  observed  changes  in  impedance  were  often  small 
making  precise  quantitative  comparison  with  theory  difficult, 
However,  certain  qualitative  conclusions  are  evident: 

1.  The  experiment  has  demonstrated  correlation  between 
measured  data  and  a  conductivity  model  which  accounts  for 
the  inertia  of  the  charge  carrier  to  within  1%  .      The 
absolute  measurement  error  (see  Table  3)  increased  with 
frequency,  with  the  largest  error  less  than  3.5%.   The 
conductivity  of  the  test  solutions  decreased  with  increased 
frequency. 

2.  The  first  measurements  of  the  non-dielectric  time 
constant,  Tc ,  defined  by  the  conductivity  model,  were 
generally  of  the  order  of  nanoseconds.   Tc  was  relatively 
insensitive  to  reagent  type  and  concentration  and  changed 
by  a  factor  of  10  for  salinity  changes  of  100. 

3.  Observed  trends  in  the  D.C.  or  real  part  of  the 
conductivity  followed  the  empirical  formula  introduced 
by  Walden  in  1906  which  characterizes  similar  data  from 
other  research. 

The  proposed  conductivity  model  has  been  shown  valid 
in  the  frequency  region  between  10  kHz  and  10  MHz,  and 
provides  a  macroscopic  rather  than  a  microscopic  description 
of  the  actual  chemical  processes  involved. 
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Future  research  is  necessary  to  provide  precise  numerical 
values  and  to  supplement  a  more  sophisticated  theoretical 
model.   This  work  directly  invites  a  number  of  follow-on 
investigations  of  phenomena  associated  with: 

-  the  temperature  dependence  of  the  conductivity,  which 
is  important  both  for  a  general  understanding  of  the 
electrolytic  conductivity  and  in  specific  areas  such  as 
biological  systems 

-  the  concentration  dependence  of  the  conductivity, 

in  the  low  concentration  regime  utilizing  more  sensitive 
test  cells  which  would  provide  more  information  on  the 
microscopic  nature  of  charge  transport  in  solutions 

-  attenuation  studies  utilizing  actual  signal  attentuation 
rather  than  measured  impedance  variation  to  further  verify 
conductivity  theory  by  independent  measurements 

-  conductivity  experiments  using  complex  solutions  formed 
by  combining  multiple  reagents  and  a  solvent  which  more 
closely  simulate  practical  electrolytic  systems  and 
provide  data  to  extend  theoretical  models  of  solutions. 

All  of  which  are  important  and  ultimately  necessary  to 
correctly  determine  the  conductivity  of  seawater. 
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-.  083 

7.933 

-  .  127 

4E  +  04 

12.517 

-.  040 

7.931 

-.06  0 

IE  +  05 

12.512 

-.002 

7.928 

.0  05 

2E  +  05 

12.507 

.040 

7.925 

.075 

4E+05 

12.499 

.121 

7.920 

.209 

IE  +  06 

12.488 

.341 

7.910 

.594 

2E  +  06 

12.478 

.783 

7.903 

i  .323 

4E  +  06 

12.465 

1.689 

7.891 

2.845 

IE  +  07 

12.549 

4.288 

8.  03i 

7.626 
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Hz 
iE  +  04 


S=.  i 
OHMS      DEG 


KBr  DATA 

S=i 
OHMS      DEG 


S  =  4 
OHMS      DEG 


10226. 000 

loIiiTooo 


-.2S6 

-Tiii" 


1078.940 

ioioT4io 


.012 

-Toio" 


279.827 
279~875 


-.003 

-TooT 


2E  +  04 
4E+04 


10229.700   -1.082   1080.570    -.080 
10207.400   -2.724   1080.460    -.290 


279.885 
279.888 


-.020 
~T073" 


1E  +  0S 
2E  +  05 


10144.100   -5.287 
9954.90  0  -10 . 14i 


1079.870    -.540    279.821 

i07iT690~""-iTo64~  "i-^Tisii 


-.137 
~~263" 


4E  +  05 
IE  +  06 


9145.100  -22.671 
7617.700  -38.12i 


1075.600   -2.602 
1069.070   -4.992 


279.493    -.667 

279~iii  -iTiio" 


2E  +  06 
4E  +  06 
iE+07 


5275. 000  -55.950 
2551.320  -73.887 


1054.530   -9.540 
1000.110  -22.143 


278.772   -2.432 
280.825   -5.877 


Hz 
iE  +  04 


S  =  10 
OHMS      DEG 


S=25 
OHMS      DEG 


S  =  49 
OHMS      DEG 


116. 035 


-.  016 


47.454 
47.445 


-.  040 
-.  023 


24.533 
24.527 


-.  0*74 
-.  043 


2E  +  04 
4E  +  04 


1161030 
116. 018 


-.00? 

-.  015 


47,435 
47.428 


-.  016 
-.  015 


24.522 
24.5.18 


,  023 
.  006 


1E  +  0S 
2E  +  05 


ii6,0i0 

115.988 


-.033 
-.  053 


47.416 
47.399 


015 
012 


24.512 
24.502 


.  013 
.  046 


4E  +  05 
iE  +  06 


1151952 
115.905 


-.095 

-.256 


47.376 
47,349 


-.  040 
-.  016 


24.488 
24 . 4^4 


.  119 

.3.15 


2E  +  06 
4E  +  06 
iE  +  07 


115.838 

115.761 


-.449 
-.862 


47.310 
47.817 


.  009 

,  177 


24.449 
24.679 


.675 
i  .839 


117.055   -2.031 


S  =  64 
HZ      OHMS      DEG 


OHMS 


S=i 


00 
DEG 


IE  +  04 

18 

750 

- 

.  094 

ii 

.864 

-.143 

2E  +  04 

18 

744 

- 

.053 

ii 

.859 

-.076 

4E+04 

18 

740 

- 

.  027 

11 

.855 

-.036 

1E  +  05 

18 

737 

- 

.002 

ii 

.852 

.001 

2E  +  05 

18 

732 

.026 

ii 

.847 

.043 

4E+05 

18 

734 

.075 

ii 

84i 

.130 

iE  +  06 

18 

713 

.198 

ii 

830 

.364 

2E  +  06 

18 

703 

.479 

ii 

82i 

.831 

4E+06 

18 

681 

1 

.  003 

ii 

805 

1  .797 

IE  +  07 

18 

884 

2 

619 

ii 

899 

4.607 
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APPENDIX  C 
CELL  n    DATA  AND  MODEL  PARAMETERS 


Hz 
IE+04 
2E+04 
4E  +  04 
1E  +  05 
2E  +  05 
4E+05 
IE  +  06 
2E+06 
4E+06 
IE  +  07 


S  =  i 
OHMS      DEG 


MaS04  DATA 
S=4 
OHMS      DEG 


S=iO 
OHMS      DEG 


22025.100 
22039.900 


.015 
.  045 


8143.700 
8152. i00 


.032 

Toi2 


3763.230 
3767.920 


,003 
,016 


22052.600 
22069.700 


.  085 
.204 


8161 .400 
8170 .200 


.062 

.  ois 


3771 .530 
37737890 


.022 
To33" 


22093. 100 

iilii .000 


.399 

7760 


8177.700 
8186.400 


.180 

.  iis 


3774.870 
377S7430 


.0^4 
.143 


22293.800 
22697. 000 


1.729 
2.970 


8205.200 
8236.300 


.686 
1.287 


3777. 430 
3781. 080 


.295 

.581 


24027.000    4.470   8344.800    2.488   3796.560    1.171 
29103.000   -1.348   9057.400    4.231   3931.130    2.608 


Hz 
iE  +  04 
2E+04 
4E  +  04 
1E+GS 
2E+05 
4E  +  05 
iE  +  06 
2E  +  06 
4E  +  06 
iE  +  07 


S=25 
OHMS      DEG 


S=49 
OHMS      DEG 


S  =  64 
OHMS      DEG 


1820.460  -.013 

1819. 930_  -100i 

1819.430  .005 

1819.080  .016 

1818.790  .039 

1318.620  .033 

1819.100  .168 

1820.230  .341 

1824.850  .683 

1867.900  1.591 


1214 
1215 
1214 
1214 
1214 
1213 
1214 
1213 
1214 
1235 


.990 
.060 
.560 
.420 
.  060 
.830 
.260 
.730 
.450 
.270 


.  040 
.050 

.  054 
.027 

.  07i 
.  100 
.  132 

.  ._  _  t_ 

.437 
1  .  047 


902 
902 
902 
902 
902 
9  02 
902 
902 
902 
916 


.440 
.600 
.470 
.520 
,400 
.350 
.720 
.480 
.850 
.970 


.  020 
.0  39 
.  038 
J_0i8 
.  053 
.  080 
,  104 
1188 
.367 
.905 


HZ 

S=i0  0 
OHMS      DEG 

iE+04 

705, 

,810 

.005 

2E  +  04 

706- 

,  110 

.  027 

4E+04 

706 

,i60 

.027 

1E  +  05 

706, 

230 

.  0i3 

2E  +  05 

706, 

,  190 

.  043 

4E  +  05 

706, 

100 

.  069 

iE  +  06 

706, 

300 

.090 

2E  +  06 

706, 

110 

.  174 

4E+06 

706, 

300 

.341 

iE  +  07 

716, 

840 

.830 
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NaCl  DATA 


Hz 

S=i 
OHMS 

DEG 

S=4 
OHMS 

DEG 

S=-i0 
OHMS 

DEG 

iE  +  04 

i603i 

.200 

.  070 

4406. 140 

.  007 

1903 

36  0 

-.015 

2E  +  04 

16032 

.300 

.  105 

4402.310 

.019 

1901 

220 

-.  002 

4E  +  Q4 

16028 

.900 

.  132 

4398.240 

.  02i 

1899 

160 

.  OOi 

1E  +  05 

16032 

.700 

.  17i 

4393.690 

.  024 

1897 

280 

.007 

2E  +  0S 

16033 

600 

.357 

4388.710 

.051 

1895 

470 

.020 

4E  +  0S 

16041 

.000 

.639 

4384.140 

.097 

1893 

820 

.045 

iE+06 

16101 

.300 

1 

.376 

4381 .230 

.  178 

1892 

770 

.075 

2E  +  06 

16257 

000 

2 

.531 

4379.220 

.350 

1892 

180 

.162 

4E  +  06 

16852 

.900 

4 

.496 

4388. 040 

.728 

1894 

660 

.330 

1E  +  07 

19927 

000 

3 

.670 

4512.610 

i  .604 

1931 

610 

.756 

Hz 

S=2C 
OHMS 

DEG 

S=49 
OHMS      DEG 

S  =  64 

OHMS 

DEG 

iE  +  04 

792 

540 

005 

424.212 

-.  026 

329. 

090 

-.  038 

2E  +  04 

792 

470 

028 

423.807 

.  001 

329. 

120 

-.  010 

4E  +  04 

792. 

060 

027 

423.383 

.  006 

329. 

085 

-.002 

iE+05 

791  , 

670 

Oil 

423. 042 

.  OOi 

329. 

098 

-.  OOi 

2E  +  0S 

791  , 

210 

042 

422.634 

.  022 

329. 

056 

.015 

4E  +  05 

790  , 

960 

067 

422.279 

.  042 

329. 

031 

.  033 

iE  +  06 

79i  . 

180 

085 

422. 043 

.  050 

329, 

061 

.  040 

2E  +  06 

79i  . 

070 

160 

42i  .737 

.  ii7 

329. 

050 

.  iOO 

4E+06 

791  . 

520 

312 

42i  .618 

.226 

329. 

148 

.  190 

iE+07 

803. 

630 

753 

427.338 

.508 

333. 

698 

.406 

HZ 
iI  +  04 


2E+04 
4E  +  04 


iE  +  05 
2E  +  oi" 


4E  +  05 
iE  +  06 


2E  +  06 
4E+06 
iE  +  07 


OHMS 


S  =  i 


00 
DEG 


323 
303 


061 
0  27 


i_  i I 


282 
272 


012 
006 


225 
225 


239 

192 


009 
027 


225 
225 


oo 


3 


174 
141 
140 
089 


036 
103 

19i" 
419 
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CELL  ti  DATA 
NaCl  PARAMETERS 


S_ 
S/rt 


Ci 

F 


Ri_ 
OHMS 


R 
OHMS 


H 


C2      R2 
F     OHMS 


10 

003 

83 

.600 

4438 

600 

11 .7000 

.  OOOi 

144 

.37  0 

i 

00 

010 

tm 

.500 

679 

000 

.6800 

.  0035 

8 

.000 

4 

Oi 

030 

1 

.220 

187 

000 

.4100 

.  0350 

.600 

10 

03 

43 

000 

.  033 

76 

500 

.2200 

.  0300 

.812 

25 

21 

65 

900 

.  019 

33 

600 

.  1120 

.  0480 

.492 

49 

81 

100 

000 

.007 

17 

750 

.  0610 

.  1300 

.231 

65 

40 

60 

000 

.  010 

14 

250 

.  0470 

.  1800 

.166 

103 

54 

120 

000 

.  009 

9 

070 

.  0385 

.2100 

.168 

s 

S/m 


MqS04  PARAMETERS 
Ci      Ri        R_        L_ 
F     OHMS      OHMS        H 


C2      R2 
F     OHMS 


26 

045 

6.  000 

2908 

600 

1  .5000 

.  OOOi 

26  0 

90  0 

1 

13 

Im 

500 

.375 

944 

600 

.2500 

.  0048 

13 

080 

3 

87 

6 

000 

.355 

358 

600 

.4010 

.  0300 

p 

020 

10 

06 

14 

000 

.  101 

164 

750 

.3250 

.  0510 

88  0 

25 

08 

40 

000 

.  030 

78 

690 

.  1850 

.2000 

262 

42 

38 

36 

000 

.053 

52 

390 

.  1300 

.3000 

184 

64 

33 

28 

000 

.  064 

36 

580 

.  1112 

.  1420 

292 

94 

55 

20 

000 

.  047 

28 

700 

.  0900 

.  1600 

25  0 

_s 

S/M 


KC1  PARAMETERS 
Ci      Ri        R         L       C2      R2 
F     OHMS      OHMS        H       F     OHMS 


10 

100 

1  .  000 

6378 

000 

5.  0000 

.  0500 

16 

.300 

i 

00 

4 

000 

.450 

734 

000 

.6600 

.  0300 

2 

.  140 

4 

00 

8 

000 

.124 

195 

060 

.3700 

.  0990 

.638 

10 

04 

25 

000 

.  080 

75 

040 

.1770 

.2500 

.265 

25 

26 

50 

000 

.  026 

31 

060 

.  0780 

.4500 

.  121 

49 

93 

55 

000 

.033 

16 

140 

.0390 

.8500 

.  084 

65 

60 

70 

000 

.023 

12 

450 

.0240  i 

.6000 

.050 

103 

97 

75 

000 

,  015 

7 

860 

.0260  i 

.  1000 

060 
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CELL  *i  DATA 
KBr  PARAMETERS 
S       _Ci      Ri~"      ~R~       _L      _C2_     R2_ 
S/m        F     OHMS      OHMS        H       F     OHMS 


iO 

500 

i  .240 

10000 

000 

.9000 

.  0020 

231 

.000 

i 

00 

200 

.890 

1064 

500 

.6000 

.  0100 

14 

.68  0 

4 

Oi 

3 

000 

.133 

278 

300 

.4250 

.  0650 

1 

43  0 

iO 

04 

ii 

000 

.  064 

115 

500 

.2450 

.  1700 

.465 

25 

i? 

2? 

000 

.  043 

47 

200 

.1090 

.3000 

.204 

4? 

6? 

53 

000 

.024 

24 

400 

.  0550 

.5500 

.105 

65 

20 

60 

000 

.015 

18 

640 

.  0465 

.8500 

.090 

102 

93 

50 

000 

.  017 

11 

780 

.0257  i 

. .2500 

.062 

CELL  #2  DATA 
NaCl  PARAMETERS 

S_      _Ci_     Ri        P_        L_  _C2  _R2 

S/m        F     OHMS      OHMS        H  '  F  OHMS 

25.21    1.000   1.310    78^.000   2.3500  .0300  2.120 

103.54   13.000    .098    224.400    .6350  .0700     .809 

MaS04  PARAMETERS 

S        Ci      Ri        P         L  C2      R2 

_S/m        F     OHMS      OHMS        H  F  OHMS 

25.08    2.000   1.160   1815.200   6.6000  .0150  3.690 

94.55    1.000    .900    703.000   2.1000  .0100  2.190 
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ANALYSIS  COMPUTER  PROGRAM 
10  OPTION  BASE  1 

ze   dim  f  v  ie  '  > z<  ie>  * t<  4  ■■  c<ie> 

30  DISP  "ENTER  DATA  FILE, LABEL 

FILE  NAMES" 
40  INPUT  P$,G$ 
50  RSSIGN*  1  TO  Pf 
60  ASSIGN*  3  TO  Q$ 
70  OISP  "WHICH  DATA  COLUMN  1-4" 

89  INPUT  B 

90  READ*  3,8  ;  L* 
180  FOR  1=1  TO  10 

110  READ**  1,1  ;  F C  I  >  ,  T  <  1  >,  T  <  2  >  ,  T 

•:  3  >  ,  T  <  4  ) 
120  ZCI>=T<B) 
130  NEXT  I 
140  ASSIGN*  1  TO  X 
150  ASSIGN*  3  TO  * 
160  DATA  6.912  6  865 , 6  633-6  534 

,6.408,6.309,6  221,6  166,6  0 

33,5.916 
170  REM    1356. . 1386, . 1437, . 1537 

6, . 1692, . 1951 ,  1913, . 1724,  1 

4412, . 12466 
ISO  FOP  1=1  TO  19 
190  READ  C(I> 
200  NEXT  I 

210  ON  KEY*  1,"HI   F"  GOSUB  280 
220  ON  KEY*  2, "MID  F B  GOSUB  370 
230  ON  KEY*  3,KL0  F"  GOSUB  469 
240  ON  KEY*  4,"Z"  GOSUB  560 
250  CLEAR  2  KEY  LABEL 
260  DISP  "SELECT  OPTION" 
270  GOTO  270 
280  CLEAR 
290  GOSUB  S90 
300  R=S9RC ABS< A9> ) 
310  L=SQR' A63(Al >  ) 
320  DISP  "R  IS  ",R 
~2Z^    DISP  "L  IS  "  •  L 
740  DISP  "ENTER  R,L" 
350  INPUT  R,L 
36^    RETURN 

370  clear  <?  gosub  S90 

380  R2=SQR(ABS<A8>>-R 

390  DISP  "P2  IS  ";R2 

400  FOR  I=S8  TO  S9 

410  C2  =  SQR<ABSCAl/<  < R~2-Z C I) ~2> X 

R2"2> >  ) 
420  DISP  "C2  IS  " ; C2 
4  30  NEXT  I 
440  DISP  "ENTER  C2" 
4  50  INPUT  C2>?  RETURN 
460  CLEAR  @  GOSUB  890 
4  70  PI  =SQR  (  ABS  <  A0  >  >-<t'R+R2> 
430  DISP  "Rl  IS  ";R1 
490  FOR  I=SS  TO  S9 
500  C1=SQRCABS<A1/(<<:R>R2)a2-Z(I 

>~2>*R1^2;>  >  ) 
510  DISP  "CI  13  u;Cl 
520  NEXT  I 
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539 

DISP  "ENTER  CI" 

549 

INPUT  CI 

559 

RETURN 

560 

CLEAR 

570 

DISP  " FORMAT; C 1, Rl ,R,L,C2,R2 
ii 

530 

DISP  CI  , Rl ,R,L,C2,R2 

599 

DISP  "ANY  CHANGES  Y/N" 

699 

INPUT  A* 

610 

IF  A*="N"  THEN  646 

629 

DISP  "INPUT  C1,R1,R,L,C2,R2" 

630 

INPUT  C1,R1,R,L,C2,R2 

640 

DISP  "VAL  TO  CRT'PTR  1/2" 

650 

INPUT  N 

66Q 

PRINTER  IS  N 

679 

PRINT  L* 

688 

PRINT  hF0RHRT-C1jR1jR;L*C2/R 

2" 

PRINT  C1,R1,R,L,C2,R2 

630 

799 

PRINT  "FORMAT, REF  Z , CALC  Z" 

719 

FOR  1=1  TO  10 

729 

w=2*pi*f«:  I  >/i  000000 

739 

C=Ca>*.  090001 

749 

G1=1'R1  «5  G2=l/R2 

759 

A1=G1A2-KW*C1>A2 

769 

A2=G2~2+<U*C2>~2 

779 

B=(R+G1/A1+G2/A2>^2 

789 

D=L-C1. A1-C2-A2 

739 

E1=<D/C)^2 

809 

E2  =  B^'<W*C>~2 

819 

E3=W^2t<D-l/<W-2tC; >-2 

320 

G  =  SQR  c.  ABS  <  (  E 1  +E2  >  s  <  B  +  E3  >  >  > 

839 

PRINT  Z(  I);  "  V;G 

849 

NEXT  I 

859 

DISP  "ANOTHER  RUN^WN" 

369 

INPUT  0$ 

870 

IF  0 *  = " Y "  THEN  570 

869 

RETURN 

839 

REM  WILL  DO  A  SIMPLE  LST  SQR 

S  FIT  TO  DATA 

390 

S1,S2,S3, S4=9 

310 

DISP  "ENTER  START  PT,END  PT 

FOR  SUM" 

329 

INPUT  88/ S3 

339 

S=S3-S8+t 

340 

FOR  I=S8  TO  S9 

359 

X=<2*PI*F< I >/ 1000000)^2 

360 

P=Z<I;-2 

370 

S1=S1+X 

330 

S2=S2+X^2 

330 

S3=S3+P 

1099  S4=S4+P*X 

1819  NEXT  I 

1020  Y1=S3-S1 

1839  B1=S1*S1/<S2*S> 

1949  Y2=S4^S2 

1950  M1=S^S1 

1060  H0=<Y1-Y2)/<<1-B1 >*M1 ) 

1070  A1=Y1-A9*M1 

1930  RETURN 

81 


APPENDIX  E 
SKIN  EFFECT  INDUCTANCE  [Ref.  27] 

To  determine  the  inductance,  the  current  distribution 
across  the  conductor  must  be  determined. 

Consider  a  cylindrical  conductor  of  radius  a,  length  I, 
across  sectional  area  A  and  conductivity  Ko  with  current 
flowing  along  the  long  axis.   The  current  is: 

.  f  >    T   Ber  (br)  +  jBei  (br) 
urJ    l      Ber  (ba)  +  jBei  (ba) 

where    R  =  £/AKo,     b  =  (2Trf  uoKo)l/2  =  V||g  g 

Ber(X)  =  1  -  &£±*    +  LllT.    .  .. 
(2!)2      (4!)2 

Bei(X)  =f-  i^l3  *  ^^-5-  ... 
4    (3!)2      (5!)2 

1 .   Case  I 

For  r  <  a,  a  =  .02  or  .01m,  fmax  =  10  Hz  and 
Ko(max)  =  10  s/m  (characteristic  of  an  electrolytic  solution) 


ba  =  2.81  x  10"5  -^fmax  Komax 
=  .281 

which  is  less  than  1,  allowing  all  the  higher  order  terms 
to  be  neglected.   Thus  we  can  say: 


82 


Ber(ba)  -  1  »  Ber(br) 
and 

Bei(ba)  =   .       =  0  (.02  maximum) 

Therefore  the  current  is: 

i(r)  =  Io  Y-z-fs   ■  Io 

which  imples  the  current  is  constant  across  the  conductor 
cross-section. 
2.   Case  II 

For  r  <  u,  u  =  .02  or  .01  m,  fmax  =10   Hz,  and 
Komax  =  10   s/m,  (characteristic  of  a  metalic  conductor) : 
ba  =  2.81  x  IO"5*  ^fmax  Komax 
=  281 
which  is  greater  than  1,  and  emphasizes  the  higher  order 
terms.   Thus,  the  current  will  be  contained  almost  entirely 
by  the  outer  surface  of  the  conductor. 

The  electrolytic  solutions  can  be  characterized 
by  the  development  in  Case  I.  The  inductance  (L)  can  be 
computed  from: 


hi    i2    =   fj     H2dV 


where  u  is  the  permeability  and  H  is  the  magnetic  field. 
For  a  cylindrical  conductor  of  radius  a,  length  h  and 
uniform  current  distribution: 


^  Li2    =   h  Uo    rir      ^ 


; 


(2^       2TTrdr   =    hT6rT 


Thus,  the  inductance  per  unit  length  is 


L  =  Uo  ,  5  x  10-a  M 

h    8tt  m 


Therefore  for  an  electrolytic  solution  in  this  configuration, 
with  h  =  10  m  =  .lm,  the  skin  effect  inductance  is: 


L  =  5  x  10"bH  =  5nH. 


The  zero  offset  adjustment  (ZOA)  is  performed  using 
mercury  for  the  short  circuit  portion.   Therefore,  the  skin 
effect  inductance  of  mercury  must  be  computed.   Since  the 
current  in  a  mercury  conductor  is  carried  in  the  outer 
surface,  as  shown  in  Case  II,  it  can  be  modeled  as  a  plane 
conductor  since  the  curvature  is  unimportant.   The  impedance 
of  the  strip  of  conductor  of  width  d,  and  unit  length  is: 


L  J  ' 


E  =  electric  field 

J  =  current  per  unit  length 


=  R  +  jo)L 


It  can  be  shown  that 

(1  t    U 

1   =  Kodb 


Ko  =  conductivity,  f  =  frequency 

,     503.3 
b  = 


^TKo" 
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Equating  the  imaginary  parts  of  both  expression 


i"L   =  ^KbTo-  u  =  27Tf 


the  impedance  becomes : 

t       1        3^162x10- "  „       ..  .    _ 

L  =  v^ua —  =  zzzz H  per  unit  length  per  unit 

wKobd        ,  \cv  -j*u 

d  yfKo  width 

For  a  =  .01m,  £  =  10"  -  107  Hz,  Ko  =  107  and  length  -  .lm, 

L  =    1.6nH     f  =  10*  Hz 
.05nH     £  =  107  Hz 

as  the  skin  effect  inductance  for  mercury. 
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